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ABSTRACT 



Aims. We analyse the effects of a first generation of fast rotating massive stars on the dynamical and chemical properties of globular 
clusters. 

Methods. We use stellar models of fast rotating massive stars, losing mass through a slow mechanical equatorial winds to produce 
material rich in H-burning products. We propose that stars with high Na and low O abundances (hereafter anomalous stars) are formed 
from matter made of slow winds of individual massive stars and of interstellar matter. The proportion of slow wind and of interstellar 
material is fixed in order to reproduce the observed Li-Na anticorrelation in NGC 6752. 

Results. In the case that globular clusters, during their lifetime, did not lose any stars, we found that to reproduce the observed ratio 
of normal to anomalous stars, a flat initial mass function (IMF) is needed, with typically a slope x = 0.55 (a Salpeter's IMF has 
x = 1.35). In the case that globular clusters suffer from an evaporation of normal stars, the IMF slope can be steeper: to have x = 1.35, 
about 96% of the normal stars would be lost. We make predictions for the distribution of stars as a function of their [O/Na] and obtain 
quite reasonable agreement with that one observed for NGC 6752. Predictions for the number fraction of stars with different values 
of helium, of the 12 C/ 13 C and 16 0/ 17 ratios are discussed, as well as the expected relations between values of [O/Na] and those of 
helium, of [C/N], of 12 C/ 13 C and of 16 0/ 17 0. Future observations might test these predictions. We also provide predictions for the 
present day mass of the clusters expressed in units of mass of the gas used to form stars, and for the way the present day mass is 
distributed between the first and second generation of stars and the stellar remnants. 
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1. Introduction 
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It has long been known that globular cluster stars present some 
striking anomalies in their content in light element^ while in 
all the Galactic globular clusters studied so far one finds "nor- 
mal" stars with detailed chemical composition similar to those 
of field stars of same metallicity (i.e., same [Fe/H]), one also ob- 
serves numerous "anomalous" main sequence and red giant stars 
that are simultaneously deficient (to various degrees) in C, O, 
and Mg, and enriched in N, Na, and A l (for recent reviews see 
iGratton et all 12004 ICharbonnelll2005l) . Additionally, the abun- 
dance of the fragile Li was found to be anticorrela ted with that of 
Na in turnofF stars in a coupl e of globular clusters (Pasquini et al.l 
l2005tlBonifacio et al.ll2007l) . 

It is clear now that these chemical peculiarities were inher- 
ited at birth by the low-mass stars we observe today in globular 
clusters, and that their root cause is H-burning through the CNO- 
cycle and the NeNa- and MgAl-chains that occurred in an early 
generation of more massive and faster evolving globular cluster 



starfl (see iPrantzos & Charbonnel l2006i hereafter IPC06L and 
references therein). In other words, compelling evidence leads 
us to believe that at least two generations of stars succeeded in 
all Galactic globular clusters during their infancy. The first one 
corresponds to the bulk of "normal" stars born with the pristine 
composition of the protocluster gas; these objects are those with 
the highest O and Mg and the lowest Na and Al abundances also 
found in their field contemporaries. The second generation con- 
tains the stars born out of material polluted to various degrees by 
the ejecta of more massive stars, and which present lower O and 
Mg and higher Na and Al abundances than their first generation 
counterparts. 

Any model attempting to explain the chemical properties of 
globular cluster stars should give an answer to the following 
questions: (1) Which type of stars produced the material en- 
riched in H-burning products? (2) What is the physical mech- 
anism responsible for selecting only material bearing the signa- 
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1 On the contrary, the content in heavy elements (i.e., Fe-group, a- 
elements) is fairly constant from star to star in any well-studied individ- 
ual Galactic globular cluster (with the notable exception of ll> Cen). 



2 As in field stars, the surface abundances of Li, C, and N also vary 
in globular cluster stars due to in situ evolutionary processes (i.e., first 
dredge-up and thermohaline instability on the red giant branch; see 
Charbonnel & Zahn 2007 and references therein). We will not address 
this point in the present paper. 
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Radiative winds 
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/ i \ 



[O/Na] = 3 
' (surface) 



Fig. 1. Schematic view of the evolution of fast rotating massive 
stars. The colours reflect the chemical composition of the various 
stellar regions and of the disc (see text for details), (top) During 
the main sequence, a slow outflowing equatorial disc forms and 
dominates matter ejection with respect to radiative winds, (mid- 
dle) At the beginning of central He-burning, the composition of 
the disc material spans the range in [O/Na] observed today in 
low-mass cluster stars. The star has already lost an important 
fraction of its initial mass, (bottom) Due to heavy mass loss, the 
star moves away from critical velocity and does not supply its 
disc anymore; radiatively-driven fast wind takes over before the 
products of He-burning reach the stellar surface. 



tures of H-processing? (3) Why does this process occur only in 
globular clusters^ 

On the basis of considerations on IMF of the potential "pol- 
luting stars", PC06 proposed massive stars (i.e., stars with an 
initial mass higher th an ~ 10 M Q ) as possible "culprits". In 
iDecressin et al.l d2007l hereafter iPaper Hi we extended this idea 
and developed stellar models sustaining our so-called "wind of 
fast rotating massive stars" (hereafter WFRMS) scenario that 



gives satisfactory answers to the previous questions, and that we 
recall below. 



2. The wind of fast rotating stars scenario - main 
guidelines 

In our WFRMS scenario, a first generation of stars forms from 
proto-clustergas already pre-enriched in heavy metals during the 
halo chemical evolution (see PC06 for a more extended discus- 
sion). They have thus the same initial composition as their field 
contemporaries; they present in particular the highest [O/Na] ra- 
tio we can observe today in any individual globular cluster. We 
assume that stars with masses between 0.1 and 120 M form 
at that stage, with an initial mass function (IMF) that has to be 
determined. 

Based on observations supporting the view that stars born 
in dense environments have faster rotational velocities than 
stars born in loose aggregates ( Keller! 20041: IStrom et ail [2005 ; 



stars born in loose aggregates (Keller 2UU4; Mromet al. 2UU0; 
iHuang & Giei2 006; Duftonetall2006), we further assume that 
the first generation massive stars are fast rotators. The evolu- 
tion of such an object is shown in Fig. \T\ The colours reflect the 
chemical composition of the various stellar regions: Green cor- 
responds to the initial chemical composition, blue and red are 
respectively for material loaded in H- and He-burning products. 
The typical [O/Na] value in the various stellar regions is indi- 
cated. In the protostars of first generation, [O/Na] is typically ~ 
0.6 (see Sect. 3). 

Fast rotating massive stars with typical time averaged veloc- 
ity of 500 km s on the main sequence easily reach the critical 
velocitjQ at the beginning of their evolution and remain near the 
critical limit during the rest of the main sequence and part of the 
central He-burning phase (see Paper I and Ekstroem et al. sub- 
mitted). As a consequence, during those phases they lose large 
amounts of material through a mechanical wincjfl which proba- 
bly leads to the formation of a slow outflowing Keplerian equa- 
torial disc. This is the kind o f process believed to occur in Be 
stars dPorter & Riviniusll2003h . 

The material ejected in the disc has two interesting charac- 
teristics: (1) due to rotational mixing that transports the products 
of the CNO cycle and of the Ne-Na and Mg-Al chains from the 
core to the stellar surface, it bears the signatures of H-burning 
and presents abundance patterns similar to the chemical anoma- 
lies observed in the second generation stars as indicated by the 
range of [O/NaQ (2) it is released into the circumstellar envi- 
ronment with a very low velocity and thus can easily be retained 
in the shallow potential well of the globular cluster. In the fol- 
lowing we shall call it the slow wind. 

The disc-star configuration lasts throughout the main se- 
quence for stars close to or at critical velocity; for the most mas- 
sive stars that are at the QF-limit (i.e., when the surface lumi- 
nosity is near the Eddington limit), it persists in the LBV phase. 



3 Indeed up to now the peculiar chemical patterns observed in glob- 
ular clusters, i.e., the O-Na and Mg-Al anticorrelations, have not been 
found in field stars. 



4 The critical velocity is the velocity at the equator such that the cen- 
trifugal acceleration balances the gravity. 

5 For example, our 60 M Q model loses 24.3 M Q in its slow wind. 

6 Early on the main sequence, the [O/Na] value in the slow wind is 
pristine (top panel). As mixing proceeds, matter with [O/Na] typical 
of H-bu rning (i.e., — 2) is ejected in the slow wind (middle panel). In 
Paper I, we showed that the yields of fast rotating massive stars (with 
an initial mass higher than ~ 20 M ) properly explain the O-Na anticor- 
relation observed today in Galactic globular clusters. We note however 
that in order to reproduce the lowest Mg values observed we had to 
postulate an increase of the 24 Mg(p, y) 25 Al nuclear reaction rate with 
respect to the published values. 
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When the star evolves away from the critical limit due to heavy 
mass loss, the radiatively-driven fast wind takes over and the 
disc is supplied by the star no longer. This happens during the 
central He-burning phase, before the He-burning products reach 
the stellar surface and contaminate the slow wind component 
(bottom panel of Fig. [TJ. From that moment on, the high-speed 
material ejected by the star and by the potential supernova will 
escape the globular cluster. We proposed that this filtering mech- 
anism, which retains in the potential well of the globular cluster 
only the slow stellar ejecta, is the physical mechanism responsi- 
ble for supplying the required H-processed material for forming 
the stars of the second generation. 

In the vicinity of massive stars, star formation may oc- 
cur, triggered for instance by the ionisation front. The large 
star-to-star spread in light elements observed today indicates 
that all the material ejected by the first generation massive 
stars did not have the time to be fully mixed before being re- 
cycled in the second generation. There is however clear ev- 
idence from the observed behaviour of both O and Li that 
the material ejected in the slo w wind is mixed with pris- 
tine interstellar matter (ISM ; see [K udryashov & Tutukov 1988; 
Denisenkov & Denisenkoyal 119891 Il990t lLanger et all Il993t 
Langer & Hoffman! 11995c iPrantzos et al.l 120071) . This suggests 
that star formation may have occurred in the vicinity of individ- 
ual massive stars from clumps in the slow wind interacting with 
pristine gas. As discussed before, the material ejected in the slow 
wind presents various degrees of enrichment in H-burning prod- 
ucts. At the beginning of the disc-star phase, pristine material is 
ejected giving birth to stars with a composition similar to that of 
the first generation. As time proceeds, the slow wind becomes 
more and more polluted in H-burning products (see Fig. 1) and 
forms stars that are more and more "anomalous". This phase 
lasts only a few million years {e.g., 4.5 Myr for a 60 M Q star). 

At this stage, the WFRMS scenario gives promising answers 
to the questions raised in Sect. 1, but it does not address the ques- 
tion of the dynamical evolution of the cluster as a whole. This is 
what we explore in the present paper, as well as other conse- 
quences regarding the He content and the carbon and oxygen 
isotopic ratios of the second generation stars. Based on quanti- 
tative theoreti cal yields of our models for fast rotating massive 
stars (iPaper J) , we first address the question o f the IM F of the 
polluter stars. We use the method developed by PC06 which re- 
lies on the observed [O/Na] distribution to derive the amount 
of polluted material with respect to the pristine composition. In 
Sect. 3 we present the general equations giving the number ra- 
tio of anomalous to normal stars, the number fraction of stars 
born with a specific chemical composition, as well as the ra- 
tio of the mass of stellar remnants to stars in the cluster today. 
The two scenarios that we have investigated regarding the early 
evolution of a Galactic globular cluster are presented in Sect. 4. 
Application to the case of NGC 6752 is discussed in Sect. 5. 
Section 6 presents our predictions for the helium content and the 
isotopic ratios 12 C/ 13 C and 16 0/ 17 respectively. Discussion and 
conclusions are the subject of Sect. 7. 



3. General equations 

Before deriving the equations that will be necessary in this paper, 
we recall that only the stars with an initial mass lower or equal to 
~ 0.8 M are still alive today. This corresponds approximately 
to the turn-off mass of ~ 12 Gyr globular clusters. We will call 
them the long-lived (hereafter LL) stars. More massive stars has 
died previously and formed dark stellar remnants. 



We note al so that, according to the nucleosynthetic predic- 
tions of IPaper fl, the polluters that are able to provide the proper 
material for the second stellar generation consist of stars with 
initial masses > 20 M Q . We assume that all the low-speed mate- 
rial ejected by those massive stars (at break-up during the main 
sequence and during the latter Or-limit phase) is used locally to 
form the second stellar generation after having been diluted to 
various degrees with pristine ISM. We do not make any assump- 
tion regarding the composition of the stellar ejecta, but rather 
use our theoretical predictions for the composition of the slow 
wincfl. 

3.1. Number ratio of normal to anomalous stars 

The first quantity we need to determine is the number ratio of 
second generation (anomalous) stars to the first generation (nor- 
mal) ones. 

If the initial total mass of first generation stars of all masses 
is M l ^ c , then the present day mass in first generation long-lived 
stars is: 



M, 



LL 



-LL> 



(1) 



ff£ is the mass fraction locked initially in first generation long- 
lived stars; it is given by integration over the IMF in the mass 
range 0.1-0.8 M . e L L is the fraction of first generation long- 
lived stars that have escaped from the cluster over its history, 
and reflects the dynamical properties of the cluster as a whole. 
For the sake of simplicity we assume that this latter parameter is 
independent of the initial stellar mass. 

Regarding the second generation long-lived stars, we must 
consider the fact that they form out of the slow wind of first 
generation massive stars that is diluted with pristine gas. The 
present-day mass locked into those stars is then: 



M, 



2G 
LL 



M^x/ sw x(l+d)x/£x(l 



■2G 



■e 2G ) 



(2) 



/sw is the mass fraction of slow winds produced by massive 
stars, which is given by theoretical stellar models, d is the pa- 
rameter reflecting global dilution with pristine ISM; it can be 
inferred from detailed analysis of the observed abundance pat- 
terns, as discussed in Sect. 5.5. f^9 and e 2 ^ are respectively the 
initial mass fraction of the second generation long-lived stars and 
the fraction of these objects that have escaped the cluster (as for 
e LL we assume that e^L i s independent of the initial stellar mass). 

Therefore, the number ratio between these two populations 
of long-lived stars is: 



'll 

JG 



LL ' 
,2G\ 



A 



x /sw x (1 + d) x 

<^L> /LL 



■2G 
LL_ 

1G 



l-e 2G 

1 LL 



1 



,1G' 
- LL 



(3) 



where (M L °) and (M*£) are the average masses of the long-lived 
stars of first and second generation respectively. The left-hand 
part of this equation can be inferred from observational stud- 
ies, following the method proposed by IPC06I and based on the 
observed [O/Na] distribution; this point is discussed further in 
Sect. 5.6. The number of escaping stars is a free parameter that 
is related to the cluster dynamics. 

The stellar IMF of both generations influences the quantities 
/ sw , , and f*°, as well as the ratio (M l L G ) / (M^) . We follow 



7 In 1PC06I the estimation of the polluter IMF for NGC 2808 was 
based on very conservative assumptions that maximised the mass of 
the H-processed ejecta released by massive stars. In particular, all the 
mass outside the He-core was assumed to have the proper composition. 
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PC06] and assume that the IMF of the globular cluster is com- 
posite. For stars of first generation more massive than 0.8 M Q , 
the IMF is approximated by a power-law of the form: 



dN 
' dM 



(4) 



with x the slope that has to be determined (x - 1.35 be- 
ing the Sal peterl 1 19551 value). On the other hand we con- 
sider that the presen t -day l og-normal distribution derived by 
Pares ce & De Marclill (|2000) reflects the IMF of long-lived stars 
(0.1-0.8 M Q ) of both first and second generations: 



In <D(M) cc A 



log(M/M c ) 
2cr 



(5) 



where Mc = 0.33 + 0.03 is the peak-mass of the log-normal 
distribution, cr = 0.34 + 0.04 a standard deviation, and A a nor- 
malisation constant. We normalise all IMF to unity: 



Jo.i 



M<S>(M)dM = 1. 



(6) 



Since we presume the same IMF for the first and second gener- 
ation long-lived stars, one has that (M™) = <M^>. 

3.2. Chemical composition of second generation stars 

The next step consists of the derivation of the equations required 
to determine the chemical abundance distribution in the matter 
out of which the second generation forms. This distribution will 
reflect directly in the star-to-star abundance variations at their 
birth. It will also reflect the distribution in the star-to-star abun- 
dance variations that we observe today in stars that have not 
undergone any change of their surface abundances by whatever 
internal mixing process, i.e. typically for stars having not yet 
evolved through their first dredge-upQ. We recall that this matter 
consists of the slow winds of individual stars of first generation 
diluted locally and to various degrees with pristine ISM. As dis- 
cussed in Sect. 2, the degree of enrichment of the slow wind by 
the H-burning products is likely to be different at the beginning 
and the end of the evolution of the massive polluter. Let us call t 
the age of the massive star polluter, then 



X, 2G ( fl ,) = (l-a,)X, sw (f) + fl,X, init 



(7) 



with X lmt and X sw being the mass fraction of the element i 
respectively in the pristine ISM and in the slow wind of the 
polluter. X 2G corresponds to the same quantity after dilution 
of the massive star ejecta with pristine material, a, being the 
local dilution parameter (see more detail in Sect. 5.5). Let us 
call m$sn(t,M) the rate of ejection of mass in the slow wind 
at time t by a star of given initial mass M. The mass ejected 
during the time dt is myw{t,M)dt. After dilution, the matter 
available to form the second generation stars from this mass 
is msw(?)df/(l - a,). Integrating over the duration of the "slow 
wind" phase, At and over the IMF, one finds an expression for 
the matter available to form the second generation stars, 



M 2G - M 1G 
m GC ~ m GC 



rtr 



m sw (t, M) 



I — a. 



dt 4>(M)dM. (8) 



8 These stars may however have suffered from some depletion of their 
lithium surface abundance. 



In this case, the present-day mass locked into these stars is then 



M 



2G 
LL 



^gcX/ll x(l 



2G 



,2G 
LL 



(9) 



By identification between Eqs. (O and © we can obtain a rela- 
tion between the global dilution parameter d and the local one 



= _L f 120 f 

/SW J20 J A, 



m S w(t,M)- 

Ar i — at 



df<D(M)dM. 



(10) 



In order to obtain the mass available to form second generation 
stars with a given abundance of the element i equal to Xr G , the 
integral over time in Eq © has to be restrained to that portion of 
the slow wind phase, during which the wind has a given abun- 
dance of the element i equal to X sw , where X sw is the abundance 
such that, after dilution, the material made of slow wind and ISM 
will have an abundance equal to Xf G . If we call Af, this duration, 
then, the mass available to form second generation stars with a 
given abundance of the element i equal to X 2G is given by 



M™(Xl G ) = M 1 £ f 12 °( f 

J20 \Ja> 



m sw (t, M)- 

At: 1 



1 



dt <X>(M)dM, 



(ID 



and the present-day mass locked into those stars is then 



M 2 L c L (Xf G ) = Ml G c (Xf G ) x fl G x (1 - el G ). (12) 

Normalised distribution is obtained by dividing Eq. (TT2l by the 
total number, Mll, of first and second generation stars given re- 
spectively by Eq. ([TJ and Eq. (fj): 



M 2G (X 2G ) 
M LL 



M 2G (X 2G )x/ 2G x(l 



„2G 
"LL 



/ll G (1 



■e 1G ) + f sw (l+d)f 2G (l 



■e 2G ) 



(13) 



Note that the above mass ratio is also equal to a number ratio 
provided the mean mass of the low mass stars in the first and 
second generation are the same, which is the hypothesis we have 
made here. 

We will use the theoretical pr edictions of our models of fast 
rotating massive stars published in Paper I. In particular, we have 
that the stars with initial masses of 20, 40, 60, and 120 M Q ex- 
pel respectively 1.7, 12.5, 24.3, and 48.0 M Q in total in their 
slow wind. In order to convolve by the polluter IMF we interpo- 
late between those values for polluter stars with different initial 



3.3. Mass of stellar remnants 

Another important quantity is the contribution of stellar rem- 
nants (i.e., white dwarfs, neutron stars or stellar black holes) to 
the total cluster mass. Indeed first and second generation stars 
more massive than the so-called long-lived stars (i.e., with initial 
masses higher than 0.8 M Q ) already died and produced dark rem- 
nantfl We call and / 2G the mass fraction to the total stellar 
mass of respectively first and second generation remnants: 



p\20 

Jo.s 



R(M)<&(M)dM. 



(14) 



9 As mentioned in Sect. 3.2, we will assume later that only long-lived 
stars formed in the second generation. However we derive here the most 
general equations. 
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We follow IPC06I for the relations between the initial mass of a 
star and the mass of its remnant: 



I 0.446 + 0.106M for M < 1OM 
R(M) = I 1.5 for 10 < M < 25M 

3 for M > 25M 



(15) 



The total mass of the remnants in the cluster today is given by: 

M rem = MgfgSl -4°) + M^/ sw (l -4m)> (16) 

with e l r 2n and e^ m the fraction of remnants lost by the cluster 
during its evolution (these quantities take into account the loss 
of both remnant progenitors and remnants). 

Therefore the mass ratio between remnants and long-lived 
stars today is obtained by combining Eq. (HJ, (0, and ( fToT l. i.e.: 



M„ 



Mi 



-rem 

el°) + f S w(l+d)f£(l 



/re^d - O + /swd + d)f r f m (l - e£j 

) 



„2G 
-rem-' 

LL 



(17) 



4. Our two scenarios for the dynamical evolution of 
the cluster 

As some of the parameters used in the previous equations are 
poorly known and are related to general properties of globular 
clusters (i.e., their IMF and their dynamical evolution) that are 
poorly constrained, we explore two simplifying scenarios. For 
the sake of simplicity, we make the following assumptions in 
both cases. We suppose that the first stellar generation contains 
stars with initial masses between 0.1 and 120 M G . Regarding 
the second generation, we assume that it consisted only of low- 
mass long-lived stars (i.e., with an initial mass between 0.1 and 
0.8 M Q ); this latter assumption will obviously minimise the con- 
straint on the polluter IMF. 

4.1. Scenario I: no loss of stars and of stellar remnants 

The first scenario we consider (left panels in Fig.f2]i assumes that 
the cluster as a whole has retained all its stars and stellar rem- 
nants. In other words, the cluster stellar component underwent 
no evaporation of stars due for instance to tidal stripping. In this 
case one has thus 



e 1G 



*2G 
-LL 



.10 
-rem 



*2G 
-rem 



o. 



(18) 



Since we assume that only low-mass stars are produced in 
the second generation, we have 



/ L 2 £ = land/ r - =0. 



2G 



(19) 



We also assume that the slope of the IMF is the same as far as the 
first and second generation long-lived stars are concerned. As a 
result of these simplifications, the general equations derived in 
Sect. 3 become 

ng = fswd+d) 

„1G 



fl 



1G 
LL 



(20) 



for the number ratio between normal and anomalous long-lived 
stars (see Eq. |3j, 



M ll( X ? G ) 

Mi I, 



flE + /sw(l+<0 



(21) 



Scenario I 
(Cluster evolution at 
constant mass 

Formation of 
globular clusters 



Scenario I 
(Mass segregation) 





Slow winds of 
massive stars 
form new stars 




Today 




Massive stars (1st gene) 
Low-mass stars (1st gene) 



Low-mass stars (2nd gene) 



Fig. 2. Schematic illustration of the scenarios discussed in 
Sect. 4. Black circles and red stars stand for first and second gen- 
eration stars respectively. In the middle panels, the dashed sym- 
bols represent fast evolving massive stars giving birth to second 
generation stars. In both cases, only long-lived low-mass stars 
form in the second generation. In Scenario I (left panels) all the 
stars remain in the cluster and the results depend on the chosen 
IMF. In Scenario II (right panels), we assume mass segregation 
and pollution occurs only at the center of the cluster; the first 
generation long-lived stars that remained in the outer part are 
assumed to be lost over the whole dynamical evolution. 



for the number fraction of stars of second generation with a given 
composition (see Eq.[T3l, and 



M„ 



1G 
rem 



M LL 



fl?+fswd+d) 



(22) 



for the fraction of the total mass in stellar remnants (see Eq.[T7Ti. 

The free parameters in Scenario I are the slope of the IMF of 
the polluters and the dilution factor. It should be noted that the 
Scenario I with no evaporation of stars is equivalent to that de- 
rived in the case where stars of both first and second generation 
leave the cluster in a similar way (i.e. where <?!9 = e?9 = e^ m = 



,2G 
"rem 



LL LL 

but are not necessarily equal to as in Eq.[T8l). 



4.2. Scenario II: mass segregation and evaporation of stars 

In Scenario II, we assume (see right panels in Fig. |2]i that the 
massive polluters of first generation were born in the center of 
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the cluster or have migrated very rapidly towards this region. 
In this case, the second generation stars are created only in the 
central region while the external part of the cluster hosts only 
first generation long-lived low-mass stars. This strong radial dis- 
tribution is assumed to stay in place until the moment where 
supernovae sweep away the residual intra-cluster gas and thus 
strongly modify the cluster potential well. As a consequence, 
stars in the external part of the cluster do not remain bound and 

are ejected into the galactic halo 

In this case, we impose a Salpeter IMF for the polluter^ 
and search for the amount of first generation long-lived stars that 
need to be lost in order to fit the observational constraints. As a 
result of our hypothesis, one has = 1, ej^ = 0, while 
remains a free parameter. The general equations become (see 
Eq.HQJandED: 



,20 
LL 
,1G 



/ sw (l +d) 
/i!l G (1-4 G l) 



Mt. 



/. 1 . G (l- e . I V i )+/ S w(l+^) 



and 



f 1G (l 



M LL flG (1 _ e lG )+/sw(1+c/) 



(23) 



(24) 



(25) 



The free parameters in Scenario II are the number fractions of 
stars lost by the clusters and the dilution factor. 



5. The case of NGC 6752 

Armed with these equations, we will first focus on the glob- 
ular cluster NG C 6752 which has a bout the same metallicity 
([Fe/H] ~ -1.56. ICarretta et al]|2007al) and initial abundance pat- 
tern as the stellar models presented in iPaper J . Our theoret- 
ical nucleosynthetic predictions can thus be directly applied. 
Additionally, NGC 6752 is the Galactic globular cluster with the 
largest set of abundance data for numerous chemical elements 
(see references in Paper I). 



5.1. Observed ratio of anomalous to normal stars in NGC 
6752 

To determine the ratio between the first generation normal and 
second generation anomalous long-lived stars in NGC 6752, we 
use the method described in PC06 and base our estimation on 
th e rO/Nal distribution function derived from the observations 
bv lCarretta et al.l d 20 07a, see Fig. [3}. Typical dispersion errors on 
O and Na abundances in the sample of 120 stars are respectively 
<To = 0.101 dex and <x Na = 0.060 dex. The total dispersion on 

the [O/Na] ratio is thus cr 0Na = ^cr Q + cr 2 ~ = 0.117, so that 

the first generation stars span a range of A [O/Na] = 2<x 0Na 
0.23 - 1 As this value is higher than the bin width of 0.2 used 
for the distribution functions displayed in Fig. [3] first generation 
stars can be found in the two bins with the highest [O/Na] ratio. 
They represent between 12% and 14% of the total number of 
long-lived stars. In the following, we round up to 85/15 the ratio 
between second and first generation long-lived stars. 



10 It should be noted that this IMF is not valid for the whole mass 
range of polluters. However our fin dings would not chan ge significantly 
with a more realistic IMF as that of Kr oupa et al.1 i 19931) 




Fig. 3. ONa distribution function in NG C 6752 based on obser - 
vations of red giant branch stars from Carrett a~et al.l ((2007a). 
(top) Distribution fu nction in a number issue taken from 
ICarretta et al.l (l2007al) : (middle) normalised distribution func- 
tion, N* is the total number of observed stars; (bottom) inverse 
cumulative distribution function for the same stars. Average 
[O/Na] error is also indicated. Dotted lines refer to stars for 
which O and Na abundances have been determined. A sample 
of stars have determinations of the Na abundance but not of O. 
In this case the abundance of O has been inferred from the Na 
abundance using the global anticorrelation between these two el- 
ements. Inclusion of this sample of stars gives the distributions 
shown by the dashed lines. 



The percentage of first generation long-lived stars we find 
in NGC 6 752 is lower than the values found in previous stud- 
ies. |PC06| find indee d 30% of normal stars in NGC 2808 while 
ICarretta et al.l d2005l) infer that less than 1/3 of the stars are nor- 
mal on the basis of CNO abundances in NGC 6397, NGC 6752, 
and 47 Tuc. This latter study does not however use very precise 
statistics as it mixes several clusters. This is discussed further in 
Sect, mi 



5.2. IMF slope for the polluter stars in Scenario I 

In Scenario I, the number fraction of second to first generation 
stars (see Eq. |20j > depends on the slope of the IMF of the pol- 
luters and on the global dilution factor d. Figure [4] displays the 
fraction of second generation stars as a function of the IMF slope 
of the polluters and for dilution factors d = and 1 . 

A Salpeter IMF with no differential evaporation leads to a 
small number of second generation stars that total up 10-17% of 
the cluster long-lived stars depending on the dilution factor (0 or 
1). In order to have 85% of second generation long-lived stars 
in the framework of Scenario I, the IMF slope must be around 
0.4 if d = (full line). This value is a lower limit as it does 
not take into account dilution of the slow wind ejecta with some 
pristine gas. For a global dilution parameter d = 1, the IMF 
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Fig. 4. Percentage of second generation stars in the total long- 
lived stars observed today as a function of the IMF slope of 
the first generation polluters. The solid line refers to long-lived 
stars born out of the slow winds of massive stars only, while 
the dashed line corresponds to the case where dilution of the 
ejecta with pristine ISM is taken into account with a global di- 
lution factor d = 1. Dotted lines show the cases where d = 1, 
ej9 = 0.66 and 0.96. The cross indicates the point corresponding 
to a Salpeter IMF and to 85% of second generation long-lived 
stars. Horizontal lines indicate the constraints coming from the 
observed [O/Na] distribution for two well-studied Galactic glob- 
ular clusters. 

slope raises up to 0.55 (see the dashed line). Estimating in the 
same way the IMF slope of the polluters for the cluster NGC 
2808, which has a higher fraction of normal first generation stars, 
results in steeper IMF (0.60 and 0.76 for respectively d = and 
1), although still quite flat with respect to the Salpeter one. Thus, 
as already discussed in IPC06L Scenario I implies relatively flat 
IMFs for the polluters. 

5.3. Fraction of stars lost by the cluster in Scenario II 

We consider now Scenario II where we use a ISalpeterl IMF for 
first generation massive stars. A dilution factor d = 1 is used. In 
order to explain the high number of second generation anoma- 
lous stars observed today, a strong loss of first generation long- 
lived stars from the cluster is required. More precisely, we derive 
e LL ~ 0-958 (see the corresponding dotted curve in Fig. [4]). This 
means that, in the present day globular cluster, only ~ 4.2% of 
the first generation long-lived stars are still present. Note that 
this is a lower limit; indeed, if some second generation stars 
also escaped in the halo, then even more first generation stars 
must have been lost by the cluster in order to reproduce the high 
number of second generation anomalous stars observed today in 
NGC 6752. It should be noted that the unpolluted pop. II stars 
found in the halo would comes largely from lower-mass clus- 
ters that dissolved after residual gas expulsion, as pointed out by 
iKroupa & BolTvl d2002l) . 



Q I I i I \ I | | | | | | | | | | | I 

0.5 1 

slope of IMF 

Fig. 5. Mass fractions of different stellar populations at birth nor- 
malised to the mass of the first generation stars, as a function 
of the slope of the polluter IMF. Full line refers to the mas- 
sive stars (20-120 M ). Dotted lines indicate the long-lived stars 
(0.1-0.8 Mq) of first and second generations (see the labels of 
the curves). The mass fraction of second generation stars corre- 
sponds to the case without dilution. 

An intermediate case between Scenario I and II is shown in 
Fig. |4] (dotted line labelled with el9 = 0.66). There we assume 
that one third of first generation long-lived stars remain in the 
cluster. In this third case, an intermediate IMF slope of 0.85 is 
needed in order to reproduce the observed ratio of anomalous to 
normal stars in NGC 6752. 

5.4. Amount of stellar residues 

Figure. [5] displays the mass fraction of different stellar popula- 
tions at birth (i.e. before any evaporation of stars) as a function 
of the IMF slope for the polluters. Normalisation in made with 
the total mass of gas used to form stars (i.e., the gas used to form 
the first stellar generation added to the one used for the dilution 
process at the birth of anomalous stars). As the IMF slope de- 
creases, massive stars initially dominate the cluster mass; for an 
IMF of ~ 0.55, the matter ejected in the slow winds represents 
~ 20% of the total mass of all first generation stars and more 
than three times the mass locked into long lived first generation 
stars. 

The left panel of Fig. [6] shows how the total stellar mass of 
the present day cluster varies as a function of the polluter IMF 
in the frame of Scenario I (continuous line). The mass is nor- 
malised to the mass of gas used to form stars. The following 
remarks can be made: 

- We see that the stellar component in the present day cluster 
represents about half of the total mass of gas used to form 
stars (more precisely between ~47 and 57% depending on 
the slope of the IMF for the polluters). The rest of the mass 
consists of material ejected by the stars under the form of fast 
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Fig. 6. Variation of the present day total mass of the cluster (con- 
tinuous line) as a function of the slope of the IMF (Scenario I, 
left panel) and as a function of the fraction of stars of the first 
generation which have been lost by the cluster, e[£ (Scenario II, 
right panel). The vertical extension of each of the three labeled 
zones gives the mass fraction of the cluster locked into long lived 
stars of first and second generation stars and in remnants. The 
line labeled as e^ n = and e l r ^ n = e!9 indicate respectively the 
case where no remnant are lost by the clusters and where they 
have been lost at the same rate as long-lived stars. The vertical 
dotted line shows the ratio of second to first generation stars is 
equal to the observed value of 85/15. 



ejecta escaping the potential well of the cluster (fast winds 
of massive stars, supernova ejecta, winds of red giants and 
asymptotic giant branch stars). 

The way the total mass varies as a function of the slope of 
the IMF does appear at first sight quite surprising. Why do 
we obtain a non-monotonic behaviour? We note that if there 
had been only a unique star formation episode, the present 
day mass cluster would be lower than the mass contained 
in the gas used to form stars. As explained above, the clus- 
ter loses fast stellar ejecta. Since their proportion increases 
when the slope of the IMF is decreased, smaller masses of 
the present day clusters are obtained with flatter IMF. In our 
case, the situation is not so straightforward, since part of the 
massive star ejecta are actually released under the form of 
slow winds which are retained by the clusters and are used 
to form the second generation stars. Therefore, the increase 
of the number of massive stars (i.e. a decrease of the slope of 
the IMF) may produce either a decrease or an increase of the 
present day mass cluster, depending on how the mass lost 
by the fast ejecta is compensated by the second generation 
stars formed from the slow ejecta and from some amount of 
ISM. From Fig. |6j we see that for an IMF slope inferior to 
~ 0.75, the mass locked in second generation stars is greater 
than the mass lost through fast ejecta, while for slopes of the 
IMF superior to 0.75, the opposite situation occurs. 
The population of remnants never exceeds those of long lived 
stars. At most, it can contain about one third of the present 
day mass in the case of a Salpeter IMF. The vertical line at an 
abscissa of 0.55 corresponds to a ratio of anomalous (second 
generation) to normal (first generation) star equal to 85/15. 



The right panel of Fig.[6]shows how the present day mass of 
the cluster varies as a function of the fraction of stars of the first 
generation which escape the cluster. Since an IMF slope of 1 .35 
is chosen in that case, we obtain for e[£ = the same results 
as in Scenario I at the corresponding IMF. As already obtained 
in Fig. 2J we see that to obtain a ratio of anomalous (second 
generation) to normal (first generation) star equal to 85/15 (see 
the vertical line with an abscissa equal to 0.958), a very large 
fraction of first generation long lived stars should have been lost 
by the cluster. For this value of e^, in the case that no remnant 
has escaped (e^ m = 0), the mass of the cluster would be dom- 
inated by the stellar residues. Only when a large fraction of the 
residues also escape the cluster (e^ = e. 1( ?), would its mass be 



dominated by the stars. In that case the present day mass of the 
cluster represents less than 10% of the mass of the gas used to 
form stars. 



5.5. Amount of dilution 
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Fig. 7. Li depletion vs [Na/Fe] value. Full c i rcles i ndicate the ob- 
servations of turnoff stars by Pasq uini et al.ld2005l) in NGC 6752. 
The theoretical predictions for the Li abundance (full and dot- 
ted lines) have been decreased by 0.3 dex to take into account 
in-situ Li depl etion during the main seque nce of low-mass stars 
(see lPC06l and lCharbonnel & Primasl2"005l for more details). Full 
lines refer to Li and Na in the ejecta of individual fast rotating 
massive stars computed in Paper I; the label indicates the cor- 
responding initial stellar mass. The dotted line corresponds to 
ejecta of a 60 M Q stars computed wi th nominal reactions rates 
(model 60rB in Decr essin etai1l2007l) . We use the minimum di- 
lution factor a m j„ =0.3. 



Let us now turn to the dilution that occurs between the slow 
ejecta of massive stars and pristine gas. If all the ejecta of all 
potential polluters were fully mixed before the formation of the 
second generation stars and no dilution with ISM occurs, Eq. (0 
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Table 1. Pristine composition and mean composition of slow 
winds after convolution by an IMF for some elements, He is in 
mass fraction. 



Element 


Initial 


IMF slope (x) 






0.4 


1.35 


[Na/Fe] 


-0.28 


0.47 


0.44 


[O/Fe] 


0.32 


0.06 


0.10 


[N/Fe] 


-0.18 


0.99 


0.94 


[C/Fe] 


-0.16 


-0.56 


-0.53 


A(Li) 


2.6 


— oo 


— OO 


He 


0.245 


0.40 


0.37 



should be replaced by: 

X f = ^7120 • ( 26 ) 

J 20 m S w®(M)dM 

with, for a given individual star msw an d %i being respectively 
the total mass ejected in its slow wind and the mean mass frac- 
tion of a given element i in its ejecta. In that case second gener- 
ation stars would all have the same composition that is given in 
Table Q] for the two extreme IMF slopes found previously (i.e., 
0.4 and 1.35Q If we take this fully mixed material and then 
dilute it in a stochastic way with pristine gas, one can form sec- 
ond generation stars with composition ranging between that of 
pristine matter (in the case of very strong dilution) and that of the 
mixed ejecta (in the case without dilution). Only a small range of 
abundance variation^ can be achieved in this way. In this case, 
the theoretical [O/Fe] ratios vary by only 0.26 dex, whereas the 
observed [O/Fe] abu ndances span more than 1 dex in NGC 6752 
dCarretta et al.ll2007al) . 

As discussed in Sect. 2, there is additional compelling evi- 
dence that the slow wind of an individual star pollutes only on a 
small scale around its stellar progenitor, and that it is diluted 
locally with pristine ISM. In order to constrain quantitatively 
this local dilu tion process we use the Li variation detected in 
NGC 6752 bv lPasquini et alj d2005l) (see Fig.|7jl. It is important 
to recall that the massive star ejecta are totally "Li-free", as this 
fragile element is destroyed in these objects. We suppose that the 
matter ejected by an individual massive stars early on the main 
sequence encounters more pristine gas and is more diluted than 
the matter ejected later. As the winds are more and more en- 
riched in Na in the course of the evolution of a massive polluter, 
we expect an anticorrelation between Li and Na: first, Li-free 
matter ejected with a low (i.e., close to initial) Na content is more 
diluted in pristine matter which is Li-rich; low-mass stars with a 
relatively high Li content and a relatively low Na abundance are 
then created. Later, Li-free matter ejected with high Na abun- 
dance is diluted with less pristine material; as a consequence, the 
newly formed stars have les s Li while their Na is h igh. The Li-Na 
anticorrelation observed by Pasq uini et alj d2005l) thus provides 
a calibration of the dilution factor a in Eq. (f7j). In this paper we 

11 Note that the ejecta composition is weighted by pristine matter re- 
leased by the polluter at the beginning of the main sequence, when 
rotational-mixing did not have the time to convey H-burning products 
from the convective core to the surface. For this reason the 60 and 
120 Mq stars eject half of their slow winds with only little abundance 
variations. This explains the very similar values given in Table[TJfor the 
two different IMF slopes. 

12 We use the two standard spectroscopic notations: [X/Fe] = 
log(X/Fe) - logCXo/Feo), and A(X) = log(X/H) + 12 with X, Fe, and H 
in number density. 



use the following expression for this parameter: 

t 

a t = 1 - (1 - anito) — (27) 
f* 

with t and f» respectively the current time and the total lifetime 
during which low-velocity winds are ejected. a m i n indicates the 
lower value for the dilution and a decreases from 1 to a m ; n when 
the polluter star evolves. 

Since Li is destroyed in massive stars, we have that 

Xlf(a t ) = a f Xjf , (28) 

where is the mass fraction of Li in second generation stars 
at birth, X™ 1 is the value derived from WMAP data and the 
standard Big B ang nucleosynthesis (A(Li)=2.65 according to 
Steigman 2006). To take into account the Li depletion (of about 
0.3 de x) occurring in low-mass stars (see lCharbonnel &~P rimas 
120051) . the values of X^ are lowered by 0.3 dex. 

In Fig.Qwe superimpose our theoretical tracks for the Li-Na 
anticorrelation obtained with a value o f a m j n equal to 0.3 on the 
observed anticorrelation determined by Pasquin i et al.1 d2005l) in 
NGC 6752. All the tracks of massive stars show a clear Li-Na 
anticorrelation with a small dispersion between those tracks (see 
the continuous lines in Fig. [7]). 

The anticorrelation follo ws the same tr e nd in our models and 
in turnoff stars observed by Pasqu ini et al.l ([2005) in NGC 6752. 
A dilution factor, a m i„, around 0.3 is needed to reproduce the 
stars with the lowest Li abundance. An overall agreement (within 
the observational errors) is obtained with this value as shown 
in Fig. [7] However it should be noted that uncertainties on the 
used nuclear reaction rates can modify this anticorrelation. To 
illustrate this point we show in Fig. [7] the wind composition of 
a 60 M Q rotating star computed with a different set of reaction 
rates (dotted line, model 60rB in paper I). In this case, the Na 
enrichment largely exceeds the observed values for the same di- 
lution factor. Thus the uncertainties about the nuclear reaction 
rates prevent firmer conclusions from being reached. 

With the above prescription for a, we can use Eq. ( fTOb to 
compute the global dilution factor d. We find that d - 1.15 for 
both scenarios. The previous findings on the mass budget ob- 
tained with d — 1 are thus still valid. 

5.6. Theoretical ONa distribution 

We compute the theoretical [O/Na] distribution function with the 
parameters we have set up previously for Scenarios I and II (see 
Eqs. [2TI and [24t. A variable dilution factor a, is used with the 
parameter a m j n = 0.3 (see Eq. 1271 ). Our results are gathered in 
Fig.[8]which displays the [O/Na] distribution function. First gen- 
eration stars representing 15% of the present day distribution are 
depicted in the hatched area in Fig. [8] 

Comparing the two upper panels of Fig. [8] we see that the 
distributions are strikingly similar. How are two such different 
scenario possible? We first note that, by construction, the number 
ratio of anomalous to normal stars is the same in the two scenar- 
ios (the parameters have been chosen in order to obtain the same 
ratio of 85/15). One can make this ratio appear in the denomi- 
nator of the Eqs. d2Tb and d24l putting in evidence /sw(l + d). 
Doing this we see that whatever the scenario considered, the dis- 
tribution is proportional to the ratio of two quantities involving 
only masses released by slow winds. A decrease of the slope of 
the IMF thus produces a similar increase of both the numerator 
and the denominator. 
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Fig. 8. Distribution function of [O/Na] obtained in Scenario I 
(top) and II (bottom) shown as thick lines. Variable dilution fac- 
tor with a m j n = 0.3 is taken into account (see Eq.|27l). The label 
"1st gene" indicates the percentage of first generation long-lived 
stars remaining in the cluster throughout its evolution. Hatched 
areas indicate the first generation or "normal" stars. We also 
superimpose the obs erved distribution histogram in NGC 6752 
dCarretta et al.ll200"7ah indicated with the dotted and dashed thin 
lines which have the same meaning as in Fig. [3] 



The peak in the distribution around [O/Na] ~ -0.1 is due 
to the stars in the mass range 20-40 M which have a longer 
evolution before reaching the break-up, so that they are more 
mixed when they release their slow winds. 

The dilution with pristine gas adds to the slow winds ma- 
terial characterised with high values of [O/Na], therefore dilu- 
tion increases the values of [O/Na]. With a a m j n parameter of 0.3 
no second generation stars can be produced with a [O/Na] ratio 
lower than -1 whereas large amounts of winds generated by the 
60 and 120 M have a [O/Na] ratio lower than -1. NGC 6752 
can possibly harbour a population of super-O-poor only detected 
with a strong Na abundance (see the dashed line in Fig. [8}. For 
those stars the O abundance is not directly measured but indi- 
rectly determined through the measure of the Na abundance and 
using the Na-0 anticorrelation. If this group is real we need to 
change the minimal amount for the mixing to a m i n =* 0.03. But 
in that case, very low abundances of Li are expected. At the mo- 
ment no Li measurements have been performed at the surface of 
non-evolved stars with such a high Na content. If feasible, such 
a measurement would repr esent an interesting ex tension of the 
anticorrelation observed by Pasq uini et alj {2005). Another pos- 
sibility to explain the s e very O-poor stars has been suggested by 
iD'Antona & Ventural d2007l) . These authors explain these stars 
as resulting from a mixing occurring in the observed star itself 
along the red giant branch. This mixing would occur only in He- 
rich stars, because in these stars mixing would not be inhibited 
by /^-gradients. 

If we compare ou r results with the distribution obtained by 
ICarretta et alj d2007al) (lower panel in Fig. [8]l we can recover 
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Fig. 9. Mass ejected by slow winds with various helium contents 
for different initial mass models. 



well stars with [O/Na] in the range between -1 and 0.7 although 
the computed distribution does not recover the exact shape of 
that observed. We however emphasise that the agreement is quite 
reasonable in view of the very small number of free parameters. 
Indeed, given a set of fast rotating massive star models and im- 
posing a given value for the ratio of anomalous to normal stars 
comparable to that observed in NGC 6752, the computed distri- 
bution depends on only one parameter which is a^n- Nature is 
probably more complicated. For example, we considered only 
one rotational velocity for our massive stars, while real stars 
probably present a distribution of velocities. Some uncertainties 
pertain to the nuclear reaction rates (see the effect in Fig. |7), 
which may affect the distribution of the composition of the slow 
winds and therefore of the second generation stars. Moreover, 
as briefly indicated above, the computed distribution depends on 
the dilution factor, a parameter which is difficult to constrain pre- 
cisely and which we modelled here in a very simple way. In view 
of all these caveats, we think that the results obtained are quite 
encouraging. In the next sections, we provide some predictions 
which might be used in the future to confirm the scenarios pre- 
sented in this paper. 

6. Helium content, 12 C/ 13 C and 16 0/ 17 isotopic 
ratios 

6.1. He ejecta from massive polluters 

Figure [9] shows the composition in helium (in mass fra ction) of 
the slow winds ejected by the massive star models of iPaper \ 
The higher the initial mass of the polluter, the higher the maxi- 
mum He value in the wind. For models of 20 and 40 M the mass 
fraction of He in the wind can reach a value of 0.30 and 0.45 re- 
spectively, while it can be as high as 0.80 in more massive mod- 
els. Such a high value in the latter case is due to the combined 
action of a more vigorous internal mixing during the evolution 
of the polluter and especially during the main sequence, and to 
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a stronger mass-loss which reveals H-processed layers. For the 
120 Mq the He mass fraction in the wind is already 0.80 at H- 
exhaustion and reaches 0.90 at the end of the LBV phase. 

6.2. He content of second generation long-lived stars 

Let us now use these predictions for the He content of the slow 
wind of massive stars, in order to derive the range in the He 
content of second generation long-lived stars. 



of fast rotating massive stars can easily produce winds with such 
a value and seem to be a good candidate to be at the origin of 
these features. However we defer the discussion of this point to 
a future paper. 



6.3. Carbon isotopic ratio 
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Fig. 10. Top panel: distribution function of He (mass fraction) of 
low-mass stars for second (white area) and first (hatched area) 
generation stars. Bottom panel: cumulative distribution function 
of He in low-mass stars. 



The number fraction of second generation long-lived stars as 
a function of their helium content is shown in Fig. [10] Only the 
case for Scenario II is shown (the results would be very similar 
in the case of Scenario I). The theoretical distribution is that ex- 
pected at the birth of the second generation starQ An extended 
tail is present toward very high values (up to 0.75) with around 
12% of the low-mass stars displaying a He content higher than 
0.4. 

It should be noted that some observational features in 

[ lobular clusters such as the mult iple sequences in a> Cen 
Bedin etalj[2004b and NGC 2808 dPiottoet al.ll2007h are ex- 
plained by_an He enrichment of the low-mass stars up to H e 
0.4 dNorrisI 120041 iPiotto etal] l2005fc Ib'Antona et alj|2005h . 
Differences in the H B morphology seem also related to high 
He value (around 0.4) dD'Antona & Caloill2004tlD'Antona etaf] 



120051; ICaloi & D'Antonal 120071: iBusso et al.ll2007l) . Our model 



13 Internal mixing processes may indeed slightly change the surface 
composition of stars during their lifetime (for instance during the red 
giant phase at the first dredge-up the l2 C/ 13 C drops to about 3.5, with 
C and Li decrease accompanied with an N increase). Thus the present 
predictions can only be compared with the observed surface composi- 
tion of present day low mass stars which are still on the Main Sequence 
or at an evolutionary stage before the first dredge-up episode. 
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Fig. 11. Number fraction of stars of the second generation com- 
posed of material with various 12 C/ 13 C ratios. The hatched area 
corresponds to first generation stars. 



We do the same exercise for 12 C and 13 C. Figure QT] shows 
the number fraction of stars as a function of their 12 C/ 13 C isotopic 
ratio at birth. Only the case with Salpeter IMF (Scenario II) is 
shown as the result is almost independent of the adopted IMF. 
We obtain a bimodal distribution with a broader peak at low 
12 C/ 13 C value. The peak around 240 is only made of first gen- 
eration star{3. We note that most of the matter ejected in the 
slow wind is near CN equilibrium, i.e., it has 12 C/ 13 C ~ 3.8. 
Dilution with pristine matter shifts the corresponding peak to 
slightly higher values. Thus many second generation stars are 
expected to be born with 12 C/ 13 C around 15-20. 

Some observations of the 12 C/ 1 3 C are available in th e liter- 
ature. For subgiants in NGC 6752, ICarretta et ail d2005l) found 
values between 3 and 11 remarkably near the values we ob- 
tain for the peak centered on low 12 C/ 13 C values. Thus the ratio 
in second generation stars is found, as expected, to be low but 
not at the equilibrium. The observed values are however slightly 
smaller than the theoretical one. This might indicate a lower di- 
lution factor. Typically, to retrieve the peak around 10, we need 
to use a dilution factor a mm of 0.1. 



14 The value of 240 is taken from galactic chemical evolution from 
Prant zos"et al.l |l996). A different value around 90 would not change 
the picture: we would still obtain a bimodal distribution with a high 
value narrow peak at 90 accompanied with a broad one at lower values. 
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Fig. 12. Number fraction of stars of the second generation com- 
posed of material with various 16 0/ 17 ratios. The hatched area 
corresponds to first generation stars. 



6.4. Oxygen isotopic ratios 

Let us first recall that, during the core H-burning phase, O is 
produced in stars with initial mass < 40 M and is slightly de- 
stroyed in more massive stars, while ls O is destroyed in all mas- 
sive stars. Since in all cases, the burning rate of ls O is faster than 
that n O, the isotopic ratic[3 16 0/ 17 decreases from 2500 (orig- 
inal value) to values between 300 (20 M Q ) and 650 (120 M ). 

Figure[T2]displays the predictions for the initial composition 
of anomalous stars obtained for Scenario II with a variable dilu- 
tion factor (a m i n = 0.3). The gap between 400 and 1000 is merely 
artificial as the lowest values are entirely due to the 20 M while 
the 40 M stars start at a ratio of 1000. So the stars between 20 
and 40 M Q will likely fill this gap. We see that only a few percent 
of stars present 15 0/ 17 ratios of the order of a few hundred at 
birth. 



6.5. Expected (anti)correlations 

Figure [13] what are the expected relations between the value 
of [O/Na] in second generation stars at birth and various other 
chemical characteristics: helium content, [C/N], 12 C/ 13 C and 
16 0/ 17 ratios. These (anti)correlations correspond to those 
which should be observed in globular clusters today provided 
only non-evolved stars are used, i.e., stars which have not yet 
undergone any change of their surface abundance due to internal 
mixing processes. Many interesting points can be seen: 



IS The nuclear rates of the reactions l7 0(p, y) and 17 0(p, a) present 
large uncertainties for temperatures higher than 100 x 10 6 K which un- 
dermine the predictions for explosive nucleosynthesis. However at tem- 
peratures reached at the center of main-sequence massive stars (around 
50 x 10 6 K) those uncertainties remain small, at similar level than that 
of 16 0(p,y). 



- The relations He- [O/Na] predicted by the different initial 
mass models present very small dispersions. This indicates 
that stars observed with a given value of [O/Na] should show 
very similar helium enrichment. 

- We see that the dispersion of the [C/N] ratio at a given value 
of [O/Na] is quite small if only the most massive stars are 
considered (M > 40 M ). Due to the 20 M Q contribution, 
the dispersion can reach values of 0.4 dex in [C/N] at [O/Na] 
equal to about -0.5. In general, the higher the initial mass, 
the higher the [C/N] ratio for a given [O/Na] ratio. This re- 
flects the fact noted above that the slow wind of lower initial 
mass stars contains mainly advanced processed material (low 
[C/N] values) because they reach the critical velocity only at 
a late time during the main-sequence phase. 

- The dispersion of the 12 C/ 13 C ratio decreases for lower val- 
ues of [O/Na]. All the stars with [O/Na] ratios below about 
-0.6 are predicted to show values of the 12 C/ 13 C ratios be- 
tween 15 and 20. 

- The isotopic ratio 16 0/ 17 presents large dispersions for a 
given [O/Na] ratio on nearly the whole range of [O/Na] val- 
ues. In the case that only very massive stars (M > 60 M Q ) 
would contribute to provide the material rich in H-burning 
products, this ratio would remain above 2000. 

It would be of course interesting to study to what extent these re- 
lations depend also on the rotational velocities considered, on the 
initial metallicity and on the mass range. It is of course the hope 
that such predictions will stimulate new observations which will 
provide further constraints on the origin of (part of) the matter 
which was used to form the anomalous stars. 

7. Recall of the main results and discussion 

In this paper we have explored some consequences of the "w inds 
of fast rotating massive stars" scenario described in Pa per j 

The two scenarios recalled above involve three free param- 
eters: the slope of the IMF of the polluters, the dilution factors 
(either the global d, or the local a,), and the number fraction of 
stars of the first generation which has been lost by the cluster 
during its lifetime. 

To constrain the values of these free parameters we used 
the following observed features: the observed number ratio of 
anomalous to normal stars and the Li-Na anticorrelation ob- 
served for turnoff stars. Given a set of stellar models, this last 
feature depends only on the local dilution facto r a t . From the 
observ ed Li-Na anticorrelation in NGC 6752 by Pasquin i et al.l 
(2005) and adopting the Eq. (f28b for a,, we obtain a value of 0.3 
for fl m i„, the minimum value of the dilution factor (see Sect. 5.5). 
As expressed by the form of a,, the mixing occurs in the vicinity 
of the mass losing star. The degree of mixing varies as a func- 
tion of the mass losing star evolution. Inserting the expression 
for a, in Eq. (TTOl i. we obtain a global dilution factor d of the or- 
der of one. It means that the total mass of second generation stars 
amounts to about twice the total mass released under the form of 
slow winds by massive stars. 

Given the above values for the dilution factors, the observed 
number ratio of anomalous to normal stars depends only on the 
two remaining free parameters. In Scenario I, the cluster loses 
no stars and no stellar remnants. In this case, the IMF is con- 
strained. In Scenario II, the slope of the IMF is set to 1.35 and 
e^L i s constrained. To produce a high number of polluted stars 
(typically 5.7 times the number of normal stars as observed in 
NGC 6752) we need either (a) in Scenario I, a flat IMF with a 
slope for massive stars around 0.55 or (b) in Scenario II, that 
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Fig. 13. Upper left: Expected anti-correlation between helium and [O/Na] for second generation stars (at birth!) born from material 
ejected through slow winds of different initial mass models. Upper right: Same as the upper left panel for the correlation between 
[C/N] and [O/Na]. Lower left: Same as the upper left panel for the correlation between 12 C/ 13 C and [O/Na]. Lower right: Same as 
the upper left panel for the correlation between 16 0/ 17 and [O/Na]. 



about 96% of the unpolluted low-mass stars has been lost by the 
cluster. 

Using the two sets of free parameters obtained in Scenarios 
I and II, we produced predictions for the following features: 
the variation of the number fraction of (non-evolved) long-lived 
stars as a function of their surface values of [O/Na], of helium, of 
12 C/ 13 C and of 16 0/ 17 0, the relations between the [O/Na] ratios 
and the helium, 12 C/ 13 C, 16 0/ 17 and [C/N] ratios. Predictions 
for the mass of the present day globular clusters expressed as 
fractions of the mass of the gas which has been used to form 
stars has also been obtained, as well as how the present day mass 



is distributed among first and second long-lived generation stars, 
and stellar remnants. 

The following results have been obtained: 

- The way the number fraction of stars varies as a function 
of the surface composition is similar in the two scenarios 
provided the adopted free parameters give the same number 
ratio of anomalous to normal stars and the same dilution fac- 
tors. 

- In both scenarios, we can produce a [O/Na] anticorrelation 
in the range -1 < [O/Na] < 0.62. Quite remarkably, this 
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range corresponds to that observed in NGC 6752 if one con- 
siders only the set of stars for which both measurements of 
oxygen and sodium are available. The possible super-O-poor 
stars, whose presence is indirectly inferred in some globular 
clusters, need a lower dilution coefficient or some additional 
mixing occurring in the star itself. 

- In both scenarios with our adopted free parameters, around 
12% of the low-mass stars are expected to display a He 
content higher than 0.4 (in mass fraction). Note that He- 
enrichment is obtained keeping the sum C+N+O equals to 
its initial value as requested by the observations. 

- A tight anti-correlation is predicted between the helium 
abundance and the value of the ratio [O/Na] at the surface of 
turnoff stars. How this anti-correlation is modified by some 
internal mixing occurring during the ascent of the red giant 
branch will be studied in a forthcoming paper. 

- More than 30% of stars are expected to have surface 12 C/ 13 C 
ratios between 10 and 30. 

- The ratios 12 C/ 13 C and [O/Na] are correlated. All the stars 
with [O/Na] ratios below about -0.6 are predicted to show 
values of the 12 C/ 13 C ratios between 15 and 20. 

- Only a few percent of stars present 16 0/ 17 ratios of the order 
of a few hundred at birth. 

- A correlation with very a large dispersion is predicted be- 
tween the ratios 16 0/ 17 and [O/Na]. In the case that only 
stars more massive than about 60 M would contribute to 
the building-up of the chemical anomalies in globular clus- 
ters, then only values above 2000 would be expected for this 
ratio. 

- A correlation between [C/N] and [O/Na] is expected, whose 
dispersion increases at lower [O/Na] values. 

- In the case of Scenario I (IMF slope of 0.55, e[£ = e LL = 
e rem = e rem = 0, fl min = 0.3), the present day mass clus- 
ter represents about half of the mass which has been used to 
form stars. The rest consists of matter ejected by stars and 
lost by the clusters. The mass locked into remnants repre- 
sents 19% of the mass of the present day cluster, the first and 
second generation stars respectively 12 and 69%. 

- In the case of Scenario II (IMF slope of 1.35, ej^ = = 
0.958, fl m in = 0.3), the present day mass cluster represents 
about 9% of the mass which has been used to form stars. The 
rest consists of matter ejected by stars and of stars and stellar 
remnants lost by the cluster. The mass locked into remnants 
represents 8% of the mass of the present day cluster, the first 
and second generation stars respectively 16 and 76%. 

Obviously our two scenarios are schematic views of the real 
formation of globular clusters. We review their respective advan- 
tages and drawbacks in Sect. 7.2. 

7. 1. Comparisons with other globular clusters 

In the present study, we have provided some constraints on the 
slope of the IMF, on the degree of evaporation of the cluster and 
on the dilution factor based on observations of the globular clus- 
ter NGC 6752. We consider to what extent using other globular 
clusters, the conclusion would have been different. 

First, let us recall the fact that the dilution factors have been 
constrained on the basis of the Li-Na anticorrelation observed 
at the surface of turnoff stars in NGC 6752. In the absence of 
similar data for other clusters, we assume that the same dilution 
factor applies in all globular clusters. 

Thus for other globular clusters the other two observed fea- 
ture which can be used to constrain the slope of the IMF of the 
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Fig. 14. Observed [O/Na] distribution function in NGC 2808 
(top), NGC 6752 (middle), and NGC 6218 (bottom). Full lines 
indicate distribution function obtained with detection or limits 
of O in stars while dashed lines indicate distribution obtained by 
a global anti-correlation relationship to obtain O abundance. 



polluters and the degree of evaporation of stars are the observed 
ratio of anomalous to normal stars and the distribution of anoma- 
lous stars with given chemical abundances. 

Th e observed ratio of anomalous to normal stars in NGC 
2808 (ICarrettaetal.ll2006l). NGC 6752 (ICarrettaet all 120073) . 
and N GC 6218 dCarretta et alj l2007bl) are respectively 30% 
dPC06h . 15% and 15%. As discussed in Sect. 5.2, a higher pro- 
portion of first generation stars induces a steeper IMF in both 
scenarios (see Fig. |4|. In the case of Scenario I, the IMF slope 
remains much flatter than a Salpeter IMF even when the frac- 
tion of first generation stars amounts to 30%. In the framework 
of Scenario II, a higher proportion of first generation stars al- 
lows the cluster to lose less first generation stars, given an IMF 
slope. However, passing from 15% to 30% for the fraction of 
first generation stars will not change a lot the degree of evapora- 
tion. From Fig.|U we can roughly estimate that it will cause the 
value of d ecrease from 0.958 to about 0.85. 

Fi gureTl4l displays the observed distribution function of the 
[O/Na] ratio in NGC 2808, NGC 6752, and NGC 6218. Despite 
that the three globular clusters have similar metallicity, they 
present distinct shapes for their [O/Na] distribution function. In 
particular, NGC 2808 presents a large population of super-O- 
poor stars. In NGC 6752 the presence of super-O-poor stars is 
only inferred from the Na abundance, and NGC 6218 has al- 
most no such stars. These differences are unlikely to be due to an 
observational bias since the globular clusters have been studied 
with the same protocol. In fact they could reflect different dy- 
namical histories among the globular clusters. It should be noted 
that the present day cluster mass of NGC 2808 is higher than that 
of NGC 6752, and NGC 6752 is heavier t han NGC 6218. This is 
in line with the finding of Carretta (2006) of indications that the 
GCs with higher masses have more pronounced abundance vari- 
ations. With a similar IMF a more massive cluster will have more 
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massive stars available to create the abundance anomalies. In ad- 
dition we can expect that more massive clusters have a deeper 
potential well and so a better ability to retain the slow winds of 
massive stars. 

From the above considerations, we can thus conclude that us- 
ing other globular clusters to constrain the free parameters of our 
two scenarios would not significantly change our conclusions 
concerning the slope of the IMF or the degree of evaporation 
of stars. However, probably other parameters, such as the ini- 
tial mass of the cluster, enter into play in the global picture. We 
defer the discussion of this point to another paper where further 
aspects of the dynamical history of the clusters will be accounted 
for. 

7.2. The pros and cons of the two scenarios 

Finally we compare the two scenarios explored in this paper. 

An advantage of Scenario I is that it allows the cluster to re- 
tain its first generation of stars and thus to suffer no or negligible 
evaporation. Looking at the m ass of young globular clusters, it 
seems a reasonable hypothesis (IQstlin e t al. 2007). Indeed these 
young globular clusters present masses in the range of the ob- 
served masses of present day old galactic globular clusters. Of 
course we are not sure that the young globular clusters we see to- 
day are formed in the same way as those which gave birth to the 
old galactic globular clusters. Most of the young globular clus- 
ters are formed during galaxy merging, while at least part of the 
old galactic globular clusters may have been formed by a rapid 
collapse a proto-halo cloud. On the other hand a disadvantage 
of Scenario I is that it imposes a very flat IMF. At the moment 
there is no direct empirical evidence that the IMF depends on 
any physical parameters such as for instanc e star formation effi - 
ciency, the metallicity or the stellar density (Ma ssev et al.ll 1995b . 
Thus this is probably at the present time a drawback of this sce- 
nario. 

The advantage of Scenario II is that it allows the use of a 
standard slope for the IMF. The disadvantage however is that it 
imposes a very serious loss of first generation stars. We propose 
a reasonable explanation for why only first generation stars are 
lost based on the mass segregation of massive stars, but we must 
admit that the requirement to lose more than 95% of this gen- 
eratio n of stars might be too demanding. iBaumgardt & Kroupal 
(2007) compute a grid of globular clusters (N-body simulations) 
to study the effect of residual-gas expulsion on the survival rate 
and final properties of star clusters. They found that some sur- 
vival clusters could lose a large part of their initial mass (up to 
85%). Models of the dynamics of globular clusters with two stel- 
lar populations are still required to check whet her the strong re- 
quire ment we need can be achieved or not (e.g jDowning & Sillsl 
120071) . 

Thus, we see that whatever scenario is realised, the pecu- 
liar abundance patterns observed in globular clusters requires 
some quite extreme assumptions: either a flat IMF or a very 
strong evaporation of stars or a combination of the two aspects. 
Moreover, as explained in more detail in paper I, we have to as- 
sume that the first generation of massive stars contained a very 
high proportion of fast rotators. Is that high proportion of fast 
rotators due to the peculiar mode of star formation in dense clus- 
ters? Are there any other observational hints supporting the ex- 
istence of these fast rotators? What would be the consequences 
of this population of massive fast rotators on the frequency of 
gamma ray bursts in very young globular clusters? We leave 
open these questions for the moment and we shall probably come 
back to them in future work. A point however which appears 



worthwhile to underline here is that massive fast rotators appear 
to be interesting objects not only to explain the abundance pat- 
terns in globular clusters, but also the high N/O ratios obse rved 
in very metal poor halo stars (Chiappi ni et al.ll2005l |2006|) and 
some C-rich Ultra Metal Poor Stars dMevnet et al.ll2006t lHirschi 
120071) . 
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